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We consider the cross section limits for light dark matter candidates (m = 0.4 to 10 GeV). We 
calculate the interaction of dark matter in the crust above underground dark matter detectors and 
find that in the intermediate cross section range, the energy loss of dark matter is sufficient to fall 
below the energy threshold of current underground experiments. This implies the existence of a 
window in the dark matter exclusion limits in the micro-barn range. 
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I. INTRODUCTION 

The evidence for the existence of dark matter (DM) 
is manifold; in particular it comes from measurements 
of galactic rotational curves, our understanding of struc- 
ture formation and from the CMB spectrum. Recent 
WMAP observations determined U^^h^ = 0.135l;oo9 and 
the baryonic fraction of matter fifcil^-'^ — 0.17 ± 0.01 [1]. 

In this paper we explore the possibility that the DM 
mass is in a range close to hadronic masses, 0.4 ^ ?ti < 10 
GeV. Masses below ~ 0.4 GeV are below the thresh- 
old of direct detection DM experiments and are there- 
fore unconstrained, with the exception of axion searches. 
The region of mass to > 10 GeV is well explored today, 
up to TeV range, from strong {a ^ 10 mb) to weak 
(cr ^ 10~^pb) cross sections, and a new generation of 
experiments reaching a ^ lO^^'^^pb is planned for the 
near future. 

The mass range to, < 10 GeV has not yet been explored 
carefully. Several dark matter underground experiments 
have sufficiently low mass threshold today: the CRESST 
[2], DAMA [3], IGEX [4], COSME [5] and ELEGANT [6] 
experiments. Except for DAMA, these experiments have 
published upper limits on the cross section assuming it is 
weak, but have not addressed the case of stronger cross 
sections,^ where the approach for extracting the cross sec- 
tion limits is substantially different, as we explain below. 
Also, recent data from an X-ray experiment XQC, [8] 
proved to be valuable in constraining low mass DM, but 
limits based on the final data have not yet been derived. 
Since XQC is a space- based experiment it is especially 
suitable for exploring the higher cross section range. In 
[9] it was shown that in the low mass range the XQC 
experiment rules out Strongly Interacting DM (SIMPs, 
[10]). Dark matter with low masses and 'intermediate' 
cross sections, several orders of magnitude smaller than 
normal hadronic cross sections, remains to be fully ana- 
lyzed and that is the focus of this work. We will abbre- 
viate DM with intermediate cross section on nucleons as 
DMIC. 



^ DAMA did publish strong cross section limits in [7], but they 
were based on a dedicated experiment which had a higher mass 
threshold m > 8 GeV. 



One of the motivations for exploring the low mass 
range is the model with baryon and antibaryonic dark 
matter, which offers a simultaneous solution for the 
baryon asymmetry and dark matter problems, intro- 
duced in [11, 12]. There, it has also been shown that 
dark matter which carries baryon number hx needs to 
be lighter than 4.5 \bx\ GeV in order for the Universe to 
have zero net baryon number. In such a scenario it can be 
natural for the DM elastic cross to be in the intermediate 
range, [11, 12]. 

Early limits from DMIC direct detection experiments 
can be found in the paper [13] by Rich, Rocchia and 
Spiro in which they reported results from a 1987 bal- 
loon experiment. Starkman et al. [14] reviewed DM con- 
straints down to a mass of 1 GeV as of 1990. Wandelt et 
al. [9] added constraints based on preliminary data from 
the more recent XQC sounding rocket experiment. The 
above constraints are discussed and updated further in 
the text. In previous works on this topic the effect of the 
halo particle's velocity distribution on the cross section 
limit was not explored. Since the only detectable light 
particles are those in the exponential tail of the velocity 
distribution, the limits on light DM are sensitive to the 
parameters in the velocity distribution, in particular to 
the value of the escape velocity cutoff. We investigate 
this sensitivity in the standard Isothermal Sphere model, 
where the DM velocity distribution function is given by 
a truncated Maxwell-Boltzmann distribution. We also 
consider the spin-independent and spin-dependent inter- 
action cases separately. Except in Section VI, we assume 
a single type of DM particle. 



II. DIRECT DARK MATTER DETECTION 

The basic principle of DM detection in underground 
experiments is to measure the nuclear recoil in elastic 
collisions, see for example [15]. The interaction of a DM 
particle of mass to. < 10 GeV, produces a recoil of a nu- 
cleus of 20 keV or less. The recoil energy (which grows 
with DM mass as in (2) below) can be measured using 
various techniques. For example, in the CRESST ex- 
periment, AI2O3 crystal was used and the temperature 
rise from the excitation of phonons due to nucleus - DM 
scattering, was measured to extract the recoil energy. 

For a given velocity distribution f{v), the differential 



rate per unit recoil energy E^ in (kg day keV) ^ in the 
detector can be expressed as 

"^^ Nt nx T" dv \v\ f{v) giv)^, (1) 
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where nx is the number density of DM particles, iV^ is 
the number of target nuclei per kg of target, ax a is the 
energy dependent scattering cross section of DM on a 
nucleus with mass number A, g{v) is the probability that 
a particle with velocity v deposits an energy above the 
threshold Eth in the detector, and Vmin is the minimum 
speed the DM particle can have and produce an energy 
deposit above the threshold. The recoil energy of the 
nucleus is given by 
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where Qcu is the scattering angle in the DM-nucleus 
center of mass frame. We will assume isotropic scattering 
as is expected at low energies. So, for instance, for A = 
16, TO = 1 GeV and an energy threshold of 600 eV, the 
minimal DM velocity to produce a detectable recoil is 
Vrain = 680 km/s, in the extreme tail of the DM velocity 
distribution. 

In order to compare cross section limits from differ- 
ent targets we will normalize them to the proton-DM 
cross section, axp- For the simplest case of interactions 
which are independent of spin and the same for protons 
and neutrons, the low energy scattering amplitude from 
a nucleus with mass number A is a coherent sum of A sin- 
gle nucleon scattering amplitudes. The matrix element 
squared therefore scales with size of nucleus as ~ ^^. In 
addition the kinematic factor in the cross section depends 
on the mass of the participants in such a way [15, 16] that 
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(3) 



where fi{A) is the reduced mass of the DM-nucleus sys- 
tem, and fj,{p) is the reduced mass for the proton-DM 
system. At higher momentum transfer q^ = 2mNER the 
scattering amplitudes no longer add in phase, and the to- 
tal cross section axAil) becomes smaller proportionally 
to the form factor F'^{q^), axA{q) — o'o-F^('Z^)- 

We take this change in the cross section into account 
when we deal with higher mass (to. > 10 GeV) dark mat- 
ter; for smaller masses the effect is negligible. We 
adopt the form factor F(g^) = exp (— l/10(qi?)^) with 
R = 1.2A^/^ fm, used also in [17, 18]. The simple expo- 
nential function is suffitiently accurate for our purposes 
and easy to implement using the Monte Carlo method to 
sample momentum transfer q, from its distribution given 
by the form factor. The procedure is described in more 
detail in Appendix B. 

For spin dependent interactions the scattering ampli- 
tude changes sign with the spin orientation. Paired nucle- 
ons therefore contribute zero to the scattering amplitude 



and only nuclei with unpaired nucleon spins are sensi- 
tive to spin dependent interactions. Due to the effect 
of coherence, the spin independent interaction is usually 
dominant, depending on the mass of the exchanged par- 
ticle [19]. Therefore, the spin dependent cross section 
limit is of interest mainly if the spin independent inter- 
action is missing, as is the case, for example, with massive 
majorana neutrinos. Another example of DM with such 
properties is photino dark matter, see [16], in the case 
when there is no mixing of left- and right- handed scalar 
quarks. The amplitude for DM-nucleus spin dependent 
interaction in the case of spin 1/2 DM, in the nonrela- 
tivistic limit, is proportional to [16, 20] 



Mr^<N\J\N> -sx 



(4) 



where J is the total angular momentum of the nucleus, 
|A^ > are nuclear states and sx is the spin of the DM 
particle. In the case of scalar DM the amplitude is 



M^<N\J\N> -{qx^) 



(5) 



where q and (f are the initial and final momenta of the 
scattering DM particle. Thus the cross section for this 
interaction is proportional to the fourth power of the ra- 
tio q/M, of DM momentum to the mass of the target 
which enters through the normalization of the wavefunc- 
tion. Therefore the spin dependent part of the interac- 
tion for scalar DM is negligible when compared to the 
spin independent part. 

We adopt the standard spin-dependent cross section 
parametrization [15] 
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where A is a parameter proportional to the spin, orbital 
and total angular momenta of the unpaired nucleon. The 
factor C is related to the quark spin content of the nu- 
cleon, C — X]^3^ij' 1 ~ u,d,s, where T^' '^ is the 
charge of the quark type q and Ag is the fraction of nu- 
cleon spin contributed by quark species q. The nuclear 
cross section normalized to the nucleon cross section is 
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The values of proton and neutron C factors, Cxp, C'xn 
vary substantially depending on the model. For targets 
of the same type - odd-n (Si, Ge) or odd-p (Al, Na, 
I) nuclei - this model dependence conveniently cancels. 
The comparison of cross sections with targets of differ- 
ent types involves the Cxp/Cxn ratio. This ratio was 
thought to have the value ^ 2 for any neutralino, based 
on the older European Muon Collaboration (EMC) mea- 
surements, but the new EMC results imply a ratio which 
is close to one for pure higgsino, and is > 10 otherwise. 
(The biggest value for the ratio is Cp/Cn ~ 500, for bino.) 
We normalize our spin dependent results to the proton 
cross section axp using Cxp/C'xn — 1 for definiteness 
below. 



In this paper we assume that the DM halo velocity dis- 
tribution is given by a truncated Maxwell-Boltzniann dis- 
tribution in the galactic reference frame, as in the Isother- 
mal Sphere Halo model [21]. We direct the z axis of the 
Earth's frame in the direction of the Local Standard of 
Rest (LSR) motion. ^ The DM velocity distribution, in 
the Earth's frame, is given by 



/(u2,wj_) = A^cxp 
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Here v^ is the local circular velocity and it is equal to 
V2 times the radial velocity dispersion in the isother- 
mal sphere model; ve is the velocity of the Earth in 
the Galactic reference frame. Throughout, superscript 
"t" indicates a tangential component. This neglects the 
Earth's motion in the radial direction which is small. 
The velocities v^ and v^ are truncated according to 
\Jv\ -V v^ < Vesc , where Vesc is the escape velocity dis- 
cussed below. 

The model above is the simplest and the most com- 
monly used model which describes a self-gravitating gas 
of coUisionless particles in thermal equilibrium. On the 
other hand numerical simulations produce galaxy halos 
which are triaxial and anisotropic and may also be locally 
clumped depending on the particular merger history (see 
[22] for a review). This indicates that the standard spher- 
ical isotropic model may not be a good approximation to 
the local halo distribution. Here we aim to extract the 
allowed DM window using the simplest halo model, but 
with attention to the sensitivity of the limit to poorly 
determined parameters of the model. The effects of the 
more detailed halo shape may be explored in a further 
work. 

We ignore here the difference between the DM veloc- 
ity distribution on the Earth, deep in the potential well 
of the solar system, and the DM velocity distribution in 
free space. This is a common assumption justified by 
Gould in [23] as a consequence of Liouville's theorem. 
Recently Edsjo et al. [24] showed that the realistic DM 
velocity distribution differs from the free space distribu- 
tion, but only for velocities f ^ 50 kni/s. Therefore, the 
free space distribution is a good approximation for our 
analysis, since for light dark matter the dominant contri- 
bution to the signal comes from high velocity part of the 
distribution. 

The velocity of the Earth in the Galactic reference 
frame is given by 



VE = VLSR + VS + VE,. 



(9) 



where vlsr is the velocity of the local standard of rest 
LSR: it moves with local circular speed in tangential di- 
rection u* QD — Vc, toward I = 90°, 6 = 0°, where I and 



LSR 



^ The Local Standard of Rest used here is the dynamical LSR, 
which is a system moving in a circular orbit around the center 
of Milky Way Galaxy at the Sun's distance. 



b are galactic longitude and latitude. The velocity of 
the Sun with respect to the LSR is vs — 16.6 km/s and 
its direction is I — 53°, b — 25° in galactic coordinates. 
VE.orb = 30 kni/s is the maximal velocity of the Earth on 
its orbit around the Sun. 

The magnitude of v'lgj^ has a considerable uncertainty. 
We adopt the conservative range Vc = (220 ± 50) km/s 
which relies on purely dynamical observations [25]. Mea- 
surements based on the proper motion of nearby stars 
give a similar central value with smaller error bars, for 
example vdRo) = (218 ± 15) km/s, from Cepheids and 
Vc{Ro) = (241 ± 17) km/s, from SgrA* (see [26] and ref- 
erences therein). The choice Vc — (220 ±50) km/s is con- 
sistent with the DAMA group analysis in [27] where they 
extracted the dependence of their cross section limits on 
the uncertainty in the Maxwellian velocity distribution 
parameters. 

Projecting the Earth's velocity on the tangential direc- 
tion (i = 90°, 5 == 0°) we get 
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where 



12 km/s; w^ = 30 cos 7 km/s where 



cos 7 — 1/2 is the cosine of the angle of the inclination 
of the plane of the ecliptic, uj — 27r/365 day~^ and to is 
June 2nd, the day in the year when the velocity of the 
Earth is the highest along the LSR direction of motion. 
In the course of the year cos[cj(i — to)] changes between 
± 1, and the orbital velocity of the Earth ranges ±15 
km/s. Taking all of the uncertainties and annual varia- 
tions into account, the tangential velocity of the Earth 
with respect to the Galactic center falls in the range 
v% = (167 to 307) km/s. 

The other parameter in the velocity distribution 
with high uncertainty is the escape velocity, Vesc = 
(450 to 650) km/s [28]. We will do our analysis with 
the standard choice of velocity distribution function pa- 
rameters, 

w|; = 230 km/s, v^ = 220 km/s, v^sc ^ 650 km/s, (11) 

and with the values of ve and Vesc from their allowed 
range, which give the lowest count in the detector and 
arc therefore most conservative: 



170 km/s, Vc = 170 km/s. 



= 450 km/s. (12) 



For experiments performed in a short time interval we 
take the value of v^ ^^^ cos[a;(t — to)] which corresponds 
to the particular date of the experiment, and the lowest 
value of v*^ allowed by the uncertainties in the value of 

Vc- 

Another effect is important in the detection of mod- 
erately interacting particles. Since particles loose energy 
in the crust rapidly (mean free path is of the order of 
100 m) only those particles which come to the detector 
from 27r solid angle above it can reach the detector with 
sufficient energy. Since the velocity distribution of the 
particles arriving to the detector from above depends on 



the detector's position on Earth with respect to the direc- 
tion of LSR motion, the detected rate for a hght DMIC 
particle will vary with the daily change of position of the 
detector. This can be a powerful signal. 



III. XQC EXPERIMENT 

For light, moderately interacting dark matter the XQC 
experiment places the most stringent constraints in the 
higher cross section range. The XQC experiment was de- 
signed for high spectral resolution observation of diffuse 
X-ray background in the 60 — 1000 eV range. The Si de- 
tector consisted of an array of 36 1 mm^ microcalorime- 
ters. Each microcalorimeter had a 7000 Angstrom layer 
of HgTe X-ray absorber. Both the HgTe and the Si lay- 
ers were sensitive to the detection. The experiment was 
performed in a 100 s flight in March, and therefore the 
Earth's velocity v% falls in the 200 to 300 km/s range. 
The experiment was sensitive to energy deposit in the 
energy range 25 — 1000 eV. For energy deposits below 25 
eV the efficiency of the detector drops off rapidly. For 
energy deposits above about 70 eV the background of X- 
rays increases, so XQC adopted the range 25-60 eV for 
extraction of DM limits, and wc will do the same. This 
translates into a conservative mass threshold for the XQC 
experiment of 0.4 GeV, obtained with Vesc = 450 km/s 
and v|; = 200 km/s, which is the lowest mass explored 
by direct DM detection apart from axion searches. 

The relationship between the number of signal events 
in the detector Ns and the scattering cross section ax a 
of DM particles on nuclei is the following 



Ns = nx f T {Nsi < vsi > asi 

+ ^Hg(< WHg > 0-Hg+ < WTo > 0-To)), 



(13) 



where iVsi and iVug are the numbers of Si and Hg (Tc) 
nuclei in the detector, nx is the local number density 
of DM particles, < wsi >, < ^Hg > and < wxe > are 
the effective mean velocities of the DM particles on the 
Si and HgTe targets, / is the efficiency of the detector, 
and T = 100 s is the data-taking time. In this energy 
range, / « 0.5. The standard value for the local DM en- 
ergy density is px — 0.3 GeV cm""^. However, numerical 
simulations combined with observational constraints [29] 
indicate that the local DM energy density px rnay have a 
lower value, 0.18 <px/(GeV cm~'^) < 0.3. In our calcula- 
tions we use both the standard value px =■ 0.3 GeV/cm'^, 
and the lower value suggested by the numerical simula- 
tions, Px = 0.2 GeV/cm'^. The cross sections ctsi, (THg, 
(Tto are calculated using equations (3) and (7). In this 
section and the next we assume that DM has dominantly 
spin-independent cross section with ordinary matter. In 
Section V we consider the case of DM which has exclu- 
sively spin-dependent cross section or when both types 
of interaction are present with comparable strength. 

XQC observed two events with energy deposit in the 
25-60 eV range, and expected a background of 1.3 events. 



The equivalent 90% cl upper limit on the number of sig- 
nal events is therefore Ng = 4.61. This is obtained by 
interpolating between 4.91 for expected background = 
1.0 event and 4.41 for expected background =1.5 events, 
for 2 observed events using table IV in ref. [30]. 

We extract the cross section limits using our simula- 
tion. Because of the screening by the Earth we consider 
only particles coming from the 27r solid angle above the 
detector, for which < fi-v >< and for them we simulate 
the interaction in the detector, for particles distributed 
according to (8). We take the direction of the LSR mo- 
tion, n as the z axis. 

We choose the nucleus i which the generated DM par- 
ticle scatters from, using the relative probability for scat- 
tering from nucleus of type i, derived in Appendix A: 



i^: = 



A 



e// 



niOxA, 



Xi J2'>T-]'^XAj 



(14) 



where Xi is the mean free path in a medium consisting of 
material with a mass number Aii Xi = {niUxAi)~^ ■ Here 
rii is the number density of target nuclei i in the crust, 
'^XAi is the scattering cross section of X on nucleus A^ 
and the effective mean free path, Ag//, is given as 



^eff = 






(15) 



In each scattering the DM particle loses energy according 
to (2), and we assume isotropic scattering in the cm. 
frame. 

We determine the effective DM velocity < wa > as 



<VA >- 



Ntot 



(16) 



where the sum is over the velocities of those DM particles 
which deposit energy in the range 25-60 eV, in a collision 
with a nucleus of type A, and Ntot is the total number of 
generated DM particles. The result depends on the angle 
between the experimental look direction, and the motion 
of the Earth. The zenith direction above the place where 
the rocket was launched, fixqci is toward h = -1-82°, I = 
75°. Thus the detector position angle compared to the 
direction of motion of the Earth through the Galaxy is 
82°. Only about 10% of the collisions have an energy 
deposit in the correct range. Putting this together gives 
the 90% confidence level XQC upper hmit on the spin 
independent cross section for DMIC shown in Figures 1 
and 2. The solid line limit is obtained using the most 
conservative set of parameters [p — 0.2 GcV/cm'^, v\, = 
200 km/s, Vesc — 450 km/s) and the dotted line is the 
limit obtained by using the standard parameter values in 
equation (11); The upper boundary of the upper domain, 
a = (10^ to 10^) mb is taken from [9]. 

When the dark matter mass is higher than 10 GeV, 
the form factor suppression of cross section is taken into 
account. We give the details of that calculation in the 
Appendix B. 



In the next section we explain how the upper bound- 
aries of the excluded region from the underground exper- 
iments shown in these figures are obtained. Also shown 
are the original limits from the balloon experiment of 
Rich et al. [13] obtained using the "standard" choices for 
DM at the time of publishing (dashed line) as well as the 
limits obtained using the conservative values of parame- 
ters (w^ = 170 km/s, since the experiment was performed 
in October, and v^sc = 450 km/s). Fig. 2 zooms in on 
the allowed window in the to < 2.4 GeV range. 



IV. UNDERGROUND DETECTION 

In this section we describe the derivation of the lower 
boundary of the DMIC window from the underground 
experiments. This is the value of axp above which the 
DM particles would essentially be stopped in the Earth's 
crust before reaching the detector. (More precisely, they 
would loose energy in the interactions in the crust and 
fall below the threshold of a given experiment.) To ex- 
tract the limit on axp we generate particles with the halo 
velocity distribution and then follow their propagation 
through the Earth's crust to the detector. We simulate 
the DM particle interactions in the detector and calculate 
the rate of the detector's response. We compare it to the 
measured rate and extract cross section limits. The basic 
input parameters of our calculation are the composition 
of the target, the depth of the detector and the energy 
threshold of the experiment. We also show the value of 
the dark matter mass threshold rriTH, calculated for the 
standard and conservative parameter values given in (11) 
and (12). The parameters are summarized in Table I for 
the relevant experiments. 



TABLE I: The parameters of the experiments used for the 
extraction of cross section limits; itith is the minimum mass 
of DM particle which can produce a detectable recoil, for 
the standard and conservative parameter choice. The energy 
threshold values E'!^^ refer to the nuclear response threshold. 
This corresponds to the electron response threshold divided 
by the quenching factor of the target. 



Experiment 


Target 


Depth 




mf}fj{mTS') 


CRESST, [2] 


AI2O3 


1400 m 


600 cV 


0.8 (1.1) GeV 


DAMA, [3] 


Nal 


1400 m 


6keV 


3.5 (5) GeV 


ELEGANT, [6] 


Nal 


442 m 


lOkeV 


5 (8) GeV 


COSME I, [5] 


Ge 


263 m 


6.5 keV 


5.5 (8) GeV 


CDMS, [31] 


Si, Go 


10.6 m 


25 keV 


9.8 (16) GeV 



where R is a uniformly distributed random number, in 
the range (0,1). After simulating the distance a parti- 
cle travels before scattering, we choose the nucleus i it 
scatters from using the relative probability as in (14). 

We take the mass density of the crust to be p — 2.7 
g/cm'^. To explore the sensitivity of the result to the 
composition of the crust we consider two different com- 
positions. First we approximate the crust as being com- 
posed of quartz, Si02, which is the most common mineral 
on the Earth and is frequently the primary mineral, with 
> 98% fraction. Then we test the sensitivity of the result 
by using the average composition of the Earth's crust: 
Oxygen 46.5 %, Silicon 28.9 %, Aluminium 8.3 % and 
Iron 4.8 %, where the percentage is the mass fraction of 
the given element. Our test computer runs showed that 
both compositions give the same result up to the first 
digit, so we used simpler composition for the computing 
time benefit. Since the DM exclusion window we obtain 
at the end of this section should be very easy to explore 
in a dedicated experiment, as we show later in the text, 
we do not aim to find precise values of the signal in the 
underground detector. 

When collisions reduce the DM velocity to less than 
the Earth's radial escape velocity, Vesc = 11 km/s, DM is 
captured by the Earth and eventually thermalized. Col- 
lisions may also reverse the DM velocity in such a way 
that the DM particle leaves the surface of the Earth with 
negative velocity: effectively, the DM particle gets re- 
flected from the Earth. The majority of light DM parti- 
cles wind up being reflected as is characteristic of diffuse 
media. The percentage of reflected particles proves not 
to depend on the cross section, as long as the mean free 
path is much smaller than the radius of the earth, but it 
does depend on DM particle mass. Light particles have 
a higher chance of scattering backward and therefore a 
higher percentage of them are reflected. The initial DM 
flux on Earth equals 2.4(1.2) 10^ (1 GeY/mx) cm'^s'^, 
taking standard (conservative) parameter values. Table 
II shows the fraction of initial flux of DM particles on 
the Earth which are captured and thermalized for vari- 
ous mass values. The fraction is, up to a percent differ- 
ence, independent of whether we make the standard or 
conservative parameter choice. 

For DM particles which are not scattered back to the 
atmosphere and which pass the depth of the detector 
before falling below the energy threshold of the detector, 
the scattering in the detector is simulated. For composite 



TABLE II: The percentage of DM particles incident on Earth 
which are captured, when Aint << Re- 



In the code, after generating particles we propagate 
those which satisfy < fi-v >< through the crust. Given 
the effective mean free path in the crust (15), the dis- 
tance traveled between two collisions in a given medium 
is simulated as 



-Ae//ln[i?] 



(17) 



mass [GeV] 


2 


4 


6 


10 


100 


thermalized [%] 


21 


30 


36 


46 


94 



targets we simulate collision with different nuclei with 
probabilities given as in (14). If the energy of the recoil 
is above Eth , we accept the event and record the velocity 



of the particle which deposited the signal. The spectral 
rate per (kg day keV) is then calculated as a sum of rates 
on the separate elements of the target, as 



dR 

dEn 



(a(t)) = E 



Ii Px < v{a{t)) >i 



A. 



i nip mx 



AE 



^XA, (18) 



where fi is the mass fraction of a given element in the 
target, px is the local DM energy density, AE is the size 
of an energy bin of a given experiment and < v{a(t)) > is 
calculated as in (16). The signal in the detector falls ex- 
ponentially with axN since the energy of DM at a given 
depth gets degraded as an exponential function of the 
cross section, see [14]. Therefore the limit on axN is 
insensitive to small changes in the rate in the detector 
coming from changes in px ; we adopt the commonly used 
value Px — 0.3 GeV cm~^ for the local DM energy den- 
sity 

We emphasize here that the spectral rate is a function 
of the relative angle a{t) between the direction of the 
motion of LSR and the position of the detector. This 
angle changes during the day as 



cos a{t) = cos (5 cos a' + sin (5 sin a' sin(ciji + (f>o) (19) 



where d is the angle between the Earth's North pole and 
the motion of LSR; a' is the angle between the position 
of the detector and the North pole, and it has a value 
of (90° - geographical latitude). The angle between the 
LSR motion and Earth's North pole is (5 = 42°, so for 
an experiment around 45° latitude (as for Gran Sasso 
experiment), a' — 45°. Therefore, in the course of a 
day, the angle between the detector and the LSR motion 
varies in the range approximately 0° to 90°. 



Fig. 3 shows the rate R per (kg 
time, (19), 



day) as a function of 



R{a{t)) - V — i^ ^^ < v{a{t)) >i axA^ (20) 

^ — ^ A; m ^ m V 



=AI,0 



Ai vfip mx 



calculated for the parameters of the CRESST experi- 
ment. We choose cio so that time t=0 corresponds to 
the moment the detector starts to move away from the 
LSR axis. We see that for these masses the rate is a 
strong function of the angle of the position of the detec- 
tor with respect to the motion of the LSR, which gives 
an interesting detection signature for detector locations 
such that this angle changes during the day. 
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Figure 3: The time dependence of the measured rate in under- 
sround detectors for m=2 GeV and m=1.5 GeV DM candidates. 



To extract our limits, we average the signal from the 
simulation dR{t)/dEfi over one day: 



< dR/dEn >-- 



1 
T 



dR{t)/dER dt. (21) 



Since the shape of the spectral rate response is a function 
of (Txp in our case (because the velocity distribution func- 
tion at the depth of detector is a function of axp due to 
the interactions in the crust) the extraction of cross sec- 
tion limits is more complicated than when the rate scales 
linearly with axp- In the region where the window is lo- 
cated, i.e. masses below 2 GeV, we perform the analysis 
based on the fit method used by the CRESST group, [2]. 
The measured spectrum is fit with an empirical function 
called B. In our case B is the sum of two falling expo- 
nentials and a constant term, since we expect the signal 
only in the few lowest energy bins. For the fit we use 
the maximum likelihood method with Poissonian statis- 
tics in each bin. The maximum likelihood of the best fit. 
Bo, is Lq. We define the background function B' as the 
difference between the best fit to the measured signal, Bo 
and some hypothesised DM signal S: B' = Bq — S. Fol- 
lowing the CRESST procedure, we set B' to zero when 
S exceeds Bq ■ When ctq is such that the simulated signal 
S is below the measured signal Bq, B' adds to the signal 
5, completing it to Bq and the likelihood is unchanged. 
With increasing tJo, when S starts to exceed the function 
Bq, B' becomes zero, and we calculate the likelihood us- 
ing S alone in such bins, leading to a new likelihood L. 
Following the CRESST collaboration prescription, (Tq ex- 
cluded at 90% CL is taken to be the value of (Tq giving 
InL — InLo = —1.28^/2 [2], since 10% of the area of a nor- 
malized Gaussian distribution of width a is 1.28(7 above 
than the peak. We show the window obtained this way 
in Figure 2 and for the low mass range, in Figure 1. 

For masses higher than 2 GeV we can use a simpler 
method, since this range of masses is already ruled out 
and our plot is only indicating from which experiments 
the constraints are coming. We calculate the response of 



the detector for different cross section values, and take 
the Hmit to be that value for which the simulated sig- 
nal is below the experiment's background. Fig 4 shows 
CRESST background together with the simulated re- 
sponse for DM particles with mass mx = 2 GcV and 
10 GeV and various values of cross section. The limits 
obtained this way for different experiments arc given in 
Figure 1. 
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Figure 4: The CRESST background and the simulated re- 
sponse of tlic detector for masses mx = 2 and mx = 10 GeV, 
and different values of spin independent cross sections uxp ■ 

The only dark matter detector sensitive to particles 
with mass < 4 GeV is CRESST. Since it is the only ex- 
periment with threshold lower than the threshold of the 
balloon experiment by Rich ct al., it extends the exist- 
ing exclusion window for intermediate cross sections. For 
the CRESST experiment we perform the calculation us- 
ing both standard and conservative parameters, because 
the size of the exclusion window is very sensitive to the 
value of mass threshold, and therefore to the parame- 
ter choice. For other underground experiments we use 
only standard parameters. In the mass range 5 < to < 10 
GeV, the ELEGANT and COSME I experiments pTace 
the most stringent constraints on a DMIC scenario, since 
they are located in shallow sites a few hundred meters 
below the ground; see Table I. Other experiments sen- 
sitive in this mass range (e.g. IGEX, COSME II) are 
located in much deeper laboratories and therefore less 
suitable for DMIC limits. We therefore present limits 
from ELEGANT and COSME I, for masses 5 to 10 GeV. 
Masses grater than 10 GeV are above the threshold of 
the CDMS experiment and this experiment places the 
most stringent lower cross section limit due to having 
the smallest amount of shielding, being only 10.6 m un- 
der ground. Fig 1 shows the cross section limits these 
experiments imply, for masses to. < 10^ GeV. 
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Figure 1: Overview of the exclusion limits for spin independent 
DM-nuclcon elastic cross section coming from the direct detection 
experiments on Earth. The limits obtained by using the conserva- 
tive parameters, as explained in the text, are plotted with a solid 
line; the dotted lines are obtained by using standard parameter val- 
ues and the dashed lines show limits published by corresponding 
experiments or in the case of XQC, by Wandelt ct al. [9]. The 
region labeled with DAM A* refers to the limits published in [7] . 
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Figure 2: The allowed window for ctJc jv ^^"^ ^ spin independent 
interaction. The region above the upper curve is excluded by the 
XQC measurements. The region below the lower curve is excluded 
by the underground CRESST experiment. The region m > 2.4 GeV 
is excluded by the experiment of Rich et al.. 



V. SPIN-DEPENDENT LIMITS 

In this section we address the case in which DM has 
a spin dependent interaction with ordinary matter. We 
consider first purely SD interaction and later we consider 
the case in which both interaction types are present. We 
focus on low masses which belong to the cross section 
window allowed by the experiment of Rich et al.. 



In this case the model dependence of the C factor ratio 
drops out in normalizing the underground experiment to 
the proton cross section. 




m^ [Gev] 

Figure 5: The allowed Spin Dependent interaction for 
(Cxp/Cxn)^ = 1- The region above the upper curve is excluded 
by XQC measurements. The region below the lower curve is ex- 
cluded by CRESST. The region m>2.4 GeV is excluded by the 
balloon experiment of Rich ct al.. 

If the DM has only a spin dependent interaction with 
ordinary matter, only the small fraction of the XQC tar- 
get with nonzero spin is sensitive to DM detection. The 
nonzero spin nuclei in the target are: Si29 (4.6 % of nat- 
ural Si), Tei25 (7 %) and Hg^gg, (16.87 %); their spin 
is due an unpaired neutron. We calculate the spin de- 
pendent cross section limits from the XQC experiment 
the same way as for the spin independent case, using the 
new composition of the target. The limiting value of the 
elastic cross section of DM with protons, <Jx^, is shown 
in Figure 5. Since the XQC target consists of n-type 
nuclei, the resulting cross section with protons is propor- 
tional to the (Cxp/Cxn)'^ factor as explained in section 
II. In Figure 4 we use the value {Cxp/Cxn)"^ = 1 which 
is the minimal value this ratio may have. We note that 
the maximal value of the ratio, based on the EMC mea- 
surements is (Cxp/Cxn)"^ = 500^ and it would shift the 
XQC limit by a factor 500^ up to higher cross sections 
(substantially extending the allowed window) . 

The spin sensitive element in the CRESST target is 
Al which has an unpaired proton in the natural isotope. 
We assume that the crust consists only of Al, since it 
is the most abundant target nucleus with non-zero spin. 




Figure 6: agj vs crgjj, for CRESST and XQC experiments, for 
mass mx = 2 GeV. The region between two curves is the allowed 
region. 

The window is extended when compared to the purely 
spin independent DM interaction, as shown in Fig. 5. 
This is mostly due to the fact that sensitive part of the 
XQC target is substantially reduced. 

In Fig. 6, for mass mx — 2 GeV, we plot the asi 
vs asD limit, assuming both types of interactions are 
present. An interesting feature in the asi vs asD de- 
pendence is that, when the spin dependent and indepen- 
dent cross sections on the target nuclei are of comparable 
magnitude, screening between two types of targets allows 
cross sections to be higher for the same rate in the de- 
tector than in the case when only one type of interaction 
is present. 



VI. CONSTRAINT ON THE FRACTION OF 
DMIC 

We now turn the argument around and use the XQC 
data to place a constraint on the fraction of allowed 
DMIC as a function of its elastic cross section. We re- 
strict consideration to values of the cross section which 
are small enough that we do not have to treat energy loss 
in the material surrounding the sensitive components of 
XQC. The maximal fraction DMIC allowed by XQC data 

cross section 



^^ can then be expressed as a function of 
using (13) 



Ns 



-,{Nsi < Wsi > crsi 



nx fT' 

+ ^Hg(< ^^Hg > 0-Hg+ < l^To > 0-To)}"^ 



(22) 



where all quantities are defined as before. 

The mass overburden of XQC can be approximated 
as [32]: A = 10"''^ g/cm , for off-angle from the cen- 
ter of the field of the detector a == (0° to 30°); A = 
10 g/cm^ for a = (30° to 100°); and A = lO'' g/cm^ 
for a > 100°. The center of the field of view points 
toward I = 90°, b = 60° which makes an angle of 32° 
with the detector position direction. Since DM par- 
ticles are arriving only from above the detector, they 
will traverse either 10 g/cm'^ or 10^^ g/cm'^ overburden. 
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Figure 7: The allowed fraction p of DM candidate as a function 
of DM-nucleon cross section. For each mass, p is calculated up to 
the values of cross sections for which the interaction in the mass 
overburden of the detector becomes important. 

For example, for values of cross section of about 0.7 
mb, m = 2 GeV DM particles start to interact in the 
10 g/cm'^ overburden, thus for cross sections above this 
value our simple approach which does not account for 
the real geometry of the detector, is not applicable any- 
more. We therefore restrict our analysis to values of the 
cross section for which neglecting the interaction in the 
overburden is a good approximation. In this domain, the 
allowed fraction of DM falls linearly with increasing cross 
section, as can be seen in equation (22) since < vdm > 
remains almost constant and is given by the halo velocity 
distribution (8). The results of the simulation are shown 
in Fig. 7, for a spin independent interaction cross sec- 
tion. An analysis valid for larger cross sections, which 
takes into account details of the geometry of the XQC 
detector, is in preparation [33]. 



periment with a detector similar to the one used in 
the XQC experiment and performed on the ground. 
Here we calculate the spectral rate of DM interac- 
tions in such detector, in order to illustrate what 
the shape and magnitude of a signal would be. 
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Figure 8: The simulated minimum rate per (kg day eV) calculated 
with axp = 2 ^b, for a DM experiment on the ground, versus de- 
posited energy En in eV, for a SI target and for a target with mass 
number A=100. The solid line indicates maximal value of the cos- 
mic ray muon background determined based on the total muon flux 
as is used in the text. 

In Fig. 8 we plot the rate per (kg-day-eV), for a Si de- 
tector and DM particle masses of mx = 1 and 2 GeV as- 
suming a spin independent interaction. In the case of an 
exclusively spin dependent interaction, the signal would 
be smaller, approximately by a factor f/A^, where /is 
the fraction of odd- nuclei in the target. The calculation 
is done for a position of a detector for which the signal 
would be the smallest. We assume a short experiment 
and do not perform averaging over time of a day because 
that would increase the signal. 

The rate scales with cross section; the rate shown in 
Fig 8 is for axp = 2 /xb, the lower limit on the cross 
section window from the XQC experiment for to = 1 
GeV. Since the unshielded muon flux on the ground is 
of the order of 2 10^ (m^ s)'^ = 2 10^ (cm^ day)"\ an 
experiment performed on the ground with an array of 
micro-calorimeter absorbers such as XQC whose target 
mass is « 100 g, should readily close this window or 
observe a striking signal. 



VII. FUTURE EXPERIMENTS 



VIII. CONCLUSION 



The window for mx '^'2A GeV in the DMIC cross 
section range could be explored in a dedicated ex- 



In this paper we have determined the limits on dark 
matter in the low mass range {m < 10 GeV) and with 
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an intermediate cross section on nucleons based on the 
final XQC data and results of underground experiments 
with low mass threshold. We also updated previous lim- 
its taking into account newer halo velocity distribution. 
We found that there is an allowed window for DM mass 
771 < 2.4 GcV and cross section a « /ib. Curiously this 
window overlaps with the mass/cross section range ex- 
pected for the H dibaryon making a possible DM candi- 
date, [12]. We showed that it should be straightforward 
experimentally to explore the window. A signal due to a 
light DMIC would have strong daily variations depending 
on the detectors position with respect to the LSR motion 
and therefore provide strong signature. 

We have benefited from discussions with D. McCam- 
mon, P. Huggins, T. Shutt, E. Sefusatti, S. Solganik, P. 
Steinhardt and L. Stodolsky. 



APPENDIX A 

The probability P{x + dx) that a particle will not scat- 
ter when propagating through a distance x + dx, equals 
the probability P{x) that it does not scatter in the dis- 
tance X, times the probability that it does not scatter 
from any type i of target nuclei in the layer dx: 

, ^ , dx 

P{x) ' 



P{x + dx) = P{x) ( 1 - XI T~ 



1 



A, 



e// 



(Al) 

By solving this differential equation one gets the proba- 
bility that a particle will travel a distance a; in a given 
medium, without scattering. 



P{x) 



^ g-2:/Ae//^ 



(A2) 



The probability for scattering once and from a given 
nuclear species i in the layer {x,x + dx), is proportional 
to the product of probabilities that a particle will not 
scatter in distance x and that it will scatter from species 
of type i in dx: 

dx 



Mx)dx - e--A=// 



A,- 



(A3) 



The probability that a particle scatters once from any 
species in a dx layer is the sum of the single particle 
probabilities "^ fi{x)dx, where 



/•OO 

/ X/,(x)dx=l. 
Jo 



(A4) 



In the simulation we want to generate the spectrum of 
distances a particle travels before scattering once from 
any of elements, using a set of uniformly distributed ran- 
dom numbers. We can achieve this by equating the dif- 
ferential probability for scattering to that of a uniformly 
distributed random number. 



fi{x) dx = dR. 



After integrating 



Z^ 



fi{x) dx 



dR, 



(A5) 



(A6) 



we get for the distribution of scattering distances x 

x = -Ae// Ini?. (A7) 

The relative frequency of scattering from a nucleus of 
type i, is then given by 



Mx)dx = ^ ^ ^2^M^, (A8) 

^i 2^ UjCrxAj 



APPENDIX B 

We assume the following function for the form factor, 
as explained in Section II, 



^2(^2)=Cxp-T^(''«)' 



(Bl) 



where q is momentum transfer and R is the nuclear ra- 
dius. For a particle moving with a given velocity v, 
the mean free path to the next collision is obtained us- 
ing the cross section atot which corresponds to (y{q) in- 
tegrated over the available momentum transfer range, 
from zero to Qmax, where Qmax = "^niNEji^rnax and 
ER.max = 2^i'^/mA{v/c)'^: 



f^tot = CTo 



j;i—FHq^)dq^ 
p- dq- 



(B2) 



After a particle travels the distance calculated from the 
mean free path described above, the collision is simu- 
lated. The momentum transfer of a collision is deter- 
mined based on the distribution given by the form factor 
function, as in the usual Monte Carlo method procedure 



^ JlPHf)^ 

j^-- F^{q^)dq^' 



(B3) 



where p is a uniformly distributed random number from 
to 1. Once the momentum transfer of the collision is 
determined, the recoil energy of the nucleus. En, and 
the scattering angle of the collision, 9cm, are uniquely 
determined. 

We repeat this procedure while following the propaga- 
tion of a particle to the detector. If the particle reaches 
the detector we simulate the collision with target nuclei. 
For each collision in the target, the energy deposited 
in the detector Er is determined as above. For each 
particle i the energy transfer determines the cross sec- 
tion with target nuclei as axAiiEu) — axA,oF^{Eii). 
The rate in the detector is found as in equation (18) 
with the only difference that in this case the sum runs 
over ^j < v(a{t))axAi >i instead of depending only on 
via{t)). 
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